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RESEARCH  IN  STOCHASTIC  PROCESSES  AND  THEIR 

APPLICATIONS 

SUMMARY  OF  RESEARCH  ACTIVITY 


Research  was  conducted  and  directed  in  the  area  of  stochastic  processes  and  their 
applications  in  engineering,  neurophysiology  and  oceanography  by  the  principal  in¬ 
vestigators,  S.  Cambanis,  G.  Kallianpur  and  M.R.  Leadbetter  and  their  associates. 
A  list  of  the  main  areas  of  research  activity  follows.  More  detailed  descriptions  of  the 
work  of  all  participants  is  given  in  the  main  body  of  the  report. 

Stochastic  differential  equations  in  infinite  dimensional  spaces 
Stochastic  differential  equation  models  for  spatially  distributed  neurons 
Propagation  of  chaos  for  interacting  systems 
Nonlinear  white  noise  analysis 
Sampling  designs  for  time  series 

Wavelets,  multiresolution  decomposition,  and  random  processes 

Non-Gaussian  stable  models:  Structure  and  inference 

Inference  for  linear  and  harmonizable  time  series 

Periodically  correlated  and  other  nonstationary  processes 

Sample  function  properties 

Random  fields  and  their  prediction 

Markov  random  field  models  for  vision 

Point  processes,  random  sets,  and  random  measures 

Random  measures  associated  with  high  levels 

Tail  inference  for  stationary  sequences 
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STAMATIS  CAMBANIS 


The  work  briefly  described  here  was  developed  in  connection  with  problems  arising 
from  and  related  to  the  statistical  communication  theory  and  the  analysis  of  stochastic 
signals  and  systems. 

Part  I  considers  questions  raised  by  the  observation  of  continuous  time  random 
signals  at  discrete  sampling  times,  and  the  transmission  or  storage  of  analog  random 
signals  in  digital  form. 

Part  II  considers  non- Gaussian  models  frequently  encountered  in  practical  appli¬ 
cations.  The  goal  is  to  learn  how  Gaussian  and  linear  signal  processing  methodologies 
should  be  adapted  to  deal  with  non- Gaussian  regimes. 

Part  III  continues  the  study  of  wavelets  and  multiresolution  analysis  for  random 
processes,  and  Part  IV  deals  with  random  filtering  and  the  harmonic  analysis  of 
nonstationary  processes. 

Item  5  is  continuing  joint  work  with  E.  Masry  of  the  University  of  California, 
San  Diego.  Items  3,  4,  6,  and  8  are  in  collaboration  with  visitors  to  the  Center 
for  Stochastic  Processes;  Houdre,  Hurd,  Fakhre-Zakeri,  Leskow,  Mandrekar,  Rosinski 
and  Surgailis.  Items  1  and  2  are  continuing  work  with  former  Ph.D.  students  Benhenni 
a-nd  Su. 


I.  DIGITAL  PROCESSING  OF  ANALOG  SIGNALS 

Continuous  time  signals  are  typically  sampled  at  discrete  times  and  inferences  are 
made  on  the  basis  of  these  samples,  which  may  be  further  quantized  (or  rounded-off) 
for  digital  processing.  Items  1  and  2  describe  work  in  progress  on  sampling  designs 
for  the  estimation  of  regression  coefficients  and  on  the  degradation  of  the  performance 
of  sampling  designs  due  to  quantization. 

1.  Sampling  designs  for  regression  coefficient  estimation  with  correlated 
errors,  [l] 

The  problem  of  estimating  regression  coefficients  from  observations  at  a  finite 
number  of  properly  designed  sampling  points  is  considered  when  the  error  process  has 
correlated  values.  Sacks  and  Ylvisaker  (1966)  found  an  asymptotically  optimal  design 
for  the  best  linear  unbiased  estimator,  which  generally  may  lack  numerical  stability 
and  requires  the  precise  knowledge  of  the  covariance  function  of  the  error  process.  Su 
and  Cambanis  (1991)  found  an  asymptotically  optimal  design  for  a  simpler  estimator 
which  is  relatively  nonparametric  (with  respect  to  the  error  covariance  function)  when 
the  error  has  no  quadratic  mean  derivative.  This  was  achieved  by  properly  adjusting 
the  median  sampling  design  and  the  simpler  estimator  introduced  by  Schoenfelder 
(1978).  Here  simpler  yet  sampling  designs  and  estimators  are  introduced  which  have 
asymptotically  optimal  performance  even  for  smoother  error  processes  (with  quadratic 
mean  derivatives). 

2.  The  effect  of  quantization  on  the  performance  of  sampling  designs.  [2] 

The  most  common  form  of  quantization  is  rounding-off,  which  occurs  in  all  digital 
systems.  A  general  quantizer  approximates  an  observed  value  by  the  nearest  among  a 
finite  number  of  representative  values.  In  estimating  weighted  integrals  of  time  series 
with  no  quadratic  mean  derivatives,  by  means  of  samples  at  discrete  times  it  is  known 
that  the  rate  of  convergence  of  the  mean  square  error  is  reduced  from  n~^  to 
when  the  samples  are  quantized  (Bucklew  and  Cambanis  (1988)).  For  smoother  time 
series,  with  k  =  1,2, .. .  quadratic  mean  derivatives,  it  is  now  shown  that  the  rate  of 
convergence  is  reduced  from  to  n~^  when  the  samples  are  quantized,  which 
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is  a  very  significant  reduction.  The  interplay  between  sampling  and  quantization  is 
also  studied,  leading  to  (asymptotically)  optimal  allocation  between  the  number  of 
samples  and  the  number  of  levels  of  quantization. 

II.  NON-GAUSSIAN  MODELS 

In  continuing  the  e.xploration  of  non-Gaussian  models  we  have  studied  a  couple 
of  stable  models.  A  new  rich  class  of  stationary  stable  processes  generalizing  moving 
averages  is  introduced  and  studied  in  Item  3,  and  the  linearity  property  of  the  pre¬ 
diction  of  heavy-tailed  autoregressive  processes  in  reversed  time  is  characterized  in 
Item  4. 


3.  Generalized  stable  moving  averages.  [3] 

No  explicit  representation  is  known  for  all  stationary  non-Gaussian  stable  pro¬ 
cesses.  The  main  two  subclasses  studied,  which  have  explicit  representations  moti¬ 
vated  by  the  Gaussian  case,  are  the  harmonizable  processes,  which  are  superpositions 
of  harmonics  with  stable  amplitudes,  and  the  moving  average  processes,  which  are 
filtered  white  stable  noise.  While  in  the  Gaussian  case,  the  latter  is  a  subclass  of 
the  former,  in  the  non-Gaussian  stable  case  the  two  classes  are  disjoint.  The  study 
of  stable  moving  average  processes  is  facilitated  by  the  fact  that  their  distribution 
is  essentially  (except  for  a  translation  and  sign)  determined  by  the  filter  function  of 
the  moving  average.  This  has  made  it  possible  to  study  distributional  properties  of 
the  process  (mixing,  ergodicity,  self-similarity,  Markov  property,  etc.)  through  the 
properties  of  the  filter  functions. 

In  this  work  the  class  of  non-Gaussian  stable  moving  average  processes  is  expanded 
substantially  by  the  introduction  of  an  appropriate  joint  randomization  of  the  filter 
function  and  of  the  stable  noise,  leading  to  stable  generalized  moving  averages  (GM.A). 
The  characterization  of  their  distribution  through  their  filter  function  and  their  mixing 
measure  leads  to  a  far  reaching  generalization  of  a  theorem  of  Ranter  (1972). 

It  is  shown  that  stable  GMA’s  contain  sums  of  independent  stable  moving  averages 
and  that  they  are  still  disjoint  from  the  harmonizable  processes,  but  are  closed  under 
time  invariant  filters,  and  that  they  are  mixing,  so  they  have  strong  ergodic  properties. 
They  lead  to  a  wealth  of  new  examples  of  self-similar  processes,  beyond  the  linear 
fractional  stable  motions,  and  also  of  processes  which  are  reflection  positive,  which  is 
a  useful  weakening  of  the  Markov  property. 

4.  The  prediction  of  heavy-tailed  autoregressive  sequences:  Regression 
versus  best  linear  prediction.  [4] 

The  prediction  of  heavy-tailed  first  order  autoregressive  sequences  is  considered. 
In  forward  time  the  regression  on  all  past  values  is  the  same  as  the  one-step  regression 
on  the  previous  value,  which  is  in  fact  linear.  In  reversed  lime  the  regression  on  all 
future  values  is  the  same  as  the  one-step  regression  on  the  immediate  future  value 
(i.e.  the  Markovian  property  is  retained)  and  we  show  that  it  is  linear  if  and  only 
if  the  innovations  have  a  semistable  distribution.  This  answers  a  question  posed  by 
Rosenblatt  (1992)  who  considered  sequences  with  finite  second  moment  and  showed 
that  regression  with  time  reversed  is  linear  if  and  only  if  the  innovations  are  Gaussian. 

When  the  distribution  of  the  innovations  is  non-Gaussian  stable,  then  both  re¬ 
gressions  in  forward  and  reversed  time  are  linear,  but  while  the  forward  regression  is 
the  best  linear  predictor,  the  regression  with  time  reversed  is  not!  The  performance 
of  linear  regression  predictor  is  compared  in  this  case  with  that  of  tlie  best  linear 
predictor. 
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III.  MULTIRESOLUTION  DECOMPOSITION  AND  WAVELET 
TRANSFORMS  OF  RANDOM  SIGNALS 

The  wavelet  approximation  of  deterministic  and  random  signals  at  given  resolution 
is  considered  in  Item  5,  which  is  a  substantial  generalization  of  earlier  work  and 
encompasses  a  much  larger  class  of  wavelets.  The  properties  of  the  wavelet  transform 
of  random  signals  are  considered  in  the  nearly  completed  work  in  Item  6.  Further 
studies  are  currently  under  way. 

5.  Wavelet  approximation  of  deterministic  and  random  signals:  Conver¬ 
gence  properties  and  rates.  [5] 

An  order  asymptotic  expansion  is  developed  for  the  error  in  the  wavelet  approx¬ 
imation  at  resolution  2“^  of  deterministic  and  of  random  signals.  The  deterministic 
signals  are  assumed  to  have  n  continuous  derivatives,  while  the  random  signals  are 
only  assumed  to  have  a  correlation  function  with  continuous  order  derivatives 
off  the  diagonal  -  a  very  mild  assumption.  For  deterministic  signals  over  the  entire 
real  line,  for  stationary  random  signals  over  finite  intervals,  and  for  nonstationary 
random  signals  with  finite  mean  energy  over  the  entire  real  line,  the  moments  of  the 
scale  function  can  be  matched  with  the  signal  smoothness  to  improve  substantially 
the  quality  of  the  approximation.  In  sharp  contrast  this  does  not  appear  to  be  gen¬ 
erally  feasible  for  nonstationary  random  signals  over  finite  intervals,  as  well  as  for 
deterministic  signals  which  are  only  locally  square  integrable. 

6.  Wavelet  transforms  of  random  processes.  [6] 

A  study  has  been  initiated  of  the  properties  of  wavelet  transforms  of  random 
processes  whose  sampled  values  appear  as  coefficients  in  the  wavelet  approximation 
of  the  process  at  a  given  resolution.  A  natural  question  is  which  properties  of  the 
process  are  inherited  to  its  wavelet  transform,  and.  conversely,  which  properties  of  the 
process  can  be  read-off  properties  of  its  wavelet  transform.  For  random  processes  with 
finite  second  moment,  properties  such  as  periodicity,  stationarity,  harmonizability, 
and  self-similarity,  are  characterized  by  means  of  analogous  properties  of  their  wavelet 
transforms  at  some  scale:  The  properties  of  the  wavelet  transform  characterize  the 
corresponding  properties  of  the  increments  of  the  process  of  order  equal  to  the  order 
of  regularity  of  the  analyzing  wavelet. 

IV.  NON-STATIONARY  PROCESSES 

In  pursuing  the  study  of  non-stationary  processes,  the  random  filters  which  pre¬ 
serve  the  normality  of  certain  non-stationary  random  inputs  are  characterized  in  Item 

7.  and  further  classes  of  non-stationary  inputs  are  currently  under  study;  earlier  work 
on  weak  laws  of  large  numbers  for  periodically  and  for  almost  periodically  correlated 
processes  which  are  not  stationary  or  harmonizable  was  substantially  revised:  and 
work  is  in  progress  jointly  with  A.G.  Miamee  of  Hampton  University  on  continuous¬ 
time  correlation-autoregressive  sequences. 

7.  Random  filters  which  preserve  the  normality  of  non-stationary  random 
inputs.  [7] 

When  a  Gaussian  signal  goes  through  a  non-random  linear  filter,  its  output  is 
also  Gaussian.  We  are  interested  in  characterising  and  identifying  those  random 
linear  filters  which  are  independently  distributed  of  their  random  inputs  and  preserve 
their  normality.  If  the  input  is  a  stationary  Gaussian  process,  then  the  output  is 
Gaussian  only  when  the  linear  filter  has  non-random  gain.  Here  we  consider  non- 
stationary  random  inputs,  for  which  the  situation  is  more  delicate.  When  the  input 
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has  stationary  independent  Gaussian  increments,  then  the  output  is  Gaussian  only 
for  linear  filters  with  either  non-random  gain  or  random  sign!  On  the  other  hand 
when  the  Gaussian  input  has  non-stationary  independent  bounded  increments,  or  is  a 
non-stationary  bounded  noise  (possibly  dependent),  or  is  harmonizable  with  diffused 
spectral  measure,  then  the  output  is  Gaussian  only  for  linear  filters  with  random  sign. 
The  non-random  characteristics  of  these  filters  can  be  identified  from  the  Gaussian 
distributions  of  the  input  and  output  processes,  and  their  random  characteristics  from 
the  joint  distribution  of  input  the  output,  which  cannot  be  Gaussian  unless  the  filter 
is  non-random. 

8.  Laws  of  large  numbers  for  periodically  and  almost  periodically  corre¬ 
lated  processes.  [8] 

This  paper  gives  results  related  to  and  including  laws  of  large  numbers  for  (possibly 
non-harmonizable)  periodically  and  almost  periodically  correlated  processes.  These 
results  admit  periodically  correlated  processes  that  are  not  continuous  in  quadratic 
mean.  The  idea  of  a  stationarizing  random  shift  is  used  to  show  that  strong  law 
results  for  weakly  stationary  processes  may  be  used  to  obtain  strong  law  results  for 
such  processes. 

.A  substantial  revision  of  this  work  from  last  year  is  being  completed.  Important 
examples  have  been  added  and  the  development  for  almost  periodically  processes  has 
been  made  simpler  and  more  transparent. 
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GOPINATH  KALLIANPUR 


As  in  recent  years,  the  major  are^is  of  my  research  have  been  the  following: 

I.  Stochastic  differential  equations  in  infinite  dimensional  spaces 

II.  Nonlinear  white  noise  analysis 

III.  Feynman  integrals  and  integration  in  Hilbert  space 

IV.  Prediction  theory  of  second  order  stationary  random  fields 

A  description  of  the  research  done  under  each  heading  is  given  below: 


I.  STOCHASTIC  DIFFERENTIAL  EQUATIONS  (SDE’s)  IN  INFINITE 
DIMENSIONAL  SPACES 

The  continuing  research  in  this  area  is  an  attempt  to  develop  a  theory  of  infi¬ 
nite  dimensional  dynamical  systems.  .4  major  emphasis  of  the  present  work  is  on 
investigating  SDE’s  in  duals  of  nuclear  spaces  driven  by  discontinuous  noise  sources, 
in  particular,  Poisson  random  measures.  Most  of  the  existing  theory  is  devoted  to 
SDE’s  or  stochastic  partial  differential  equations  (SPDE)  driven  by  cylindrical  Brow¬ 
nian  motions  or  space-time  Wiener  processes  primarily  because  of  its  mathematical 
elegance  and  the  link  with  diffusion  processes. 

However,  in  problems  of  neuronal  behavior,  environmental  pollution  and  fluid 
mechanics  (to  name  only  a  few  fields  of  application)  it  seems  more  natural  to  consider 
dynamic  models,  i.e.  SDE’s  or  SPDE’s  driven  by  Poisson  random  measures.  The 
diffusion  approximations  that  can  be  derived  from  them  throw  additional  light  on  the 
continuous  models.  The  new  work  (jointly  with  J.  Xiong)  described  below  pertains 
to  recent  work  on  the  application  to  environmental  pollution  [3],  uses  the  techniques 
and  results  of  the  following  papers  partially  described  in  the  .Annual  Scientific  Report 
for  1990-1991  and  is  now  completed: 

1.  The  existence  and  uniqueness  of  the  solution  of  nuclear  space-valued 
stochastic  differential  equations  driven  by  Poisson  random  measures  (with 
G.  Hardy,  S.  Ramasubramanian  and  J.  Xiong)  [1] 

In  this  paper,  we  study  SDE’s  on  duals  of  nuclear  spaces  driven  by  Poisson  random 
measures.  The  existence  of  a  weak  solution  is  obtained  by  the  Galerkin  method. 
For  uniqueness,  a  class  of  f^-valued  processes  which  are  called  Good  processes  is 
introduced.  An  equivalence  relation  is  established  between  SDE’s  driven  by  Poisson 
random  measures  and  those  by  Good  processes.  The  uniqueness  is  established  by 
extending  the  Yamada-Watanabe  argument  to  the  SDE's  driven  by  Good  processes. 
This  is  an  extension  to  discontinuous  infinite  dimensional  SDE's  of  work  done  by  G. 
Kallianpur,  1.  Mitoma  and  R.  Wtflpert  for  nuclear  space  valued  diffusions  [Stochastics. 
29.  1-45,  (1990)]. 

2,  Asymptotic  behavior  of  a  system  of  interacting  stochastic  differential 
equations  driven  by  Poisson  random  measures  (with  J.  -Xiong)  [2] 

In  this  paper,  we  study  a  system  of  interacting  stochastic  differential  equations 
taking  values  in  nuclear  spaces  and  driven  by  Poisson  random  measures.  We  also  con¬ 
sider  the  .McKean- Vlasov  equation  associated  with  the  .system.  We  show  that  under 
suitable  conditions  the  system  has  a  unique  solution  and  the  sequence  of  its  empirical 
distributions  converges  to  the  solution  of  the  .McKean-Vlasov  equation  when  the  size 
of  the  system  tends  to  infinity.  The  results  are  applied  to  the  voltage  potentials  of 


a  large  system  of  neurons  and  a  law  of  large  numbers  for  the  empirical  measure  is 
obtained. 

3.  Stochastic  models  of  environmental  pollution  (with  J.  Xiong)  [3] 

In  this  paper,  we  consider  several  stochastic  models  arising  from  environmental 
problems.  First,  we  study  the  pollution  in  a  domain  where  undesired  chemicals  are 
deposited  at  random  times  and  locations  according  to  Poisson  streams.  Incorpo¬ 
rated  with  drift  and  dispersion,  the  chemical  concentration  can  be  modeled  by  linear 
stochastic  partial  differential  equations  (SPDE)  which  are  solved  by  applying  the 
general  results  of  Kallianpur  and  Xiong  (SDE’s  in  infinite  dimensions:  .4  brief  survey 
and  some  new  directions.  Center  for  Stochastic  Processes  Technical  Report  No.  372. 
Sept.  92). 

We  examine  in  a  somewhat  more  general  context,  the  stochastic  dynamic  model 
considered  by  Kwakernaak  and  by  Curtain,  and  look  at  the  problem  in  the  framework 
of  general  SPDE's;  Let 


X  e  <V  =  [0,  f']'*.  X  =  (xi, - Xi). 


The  underlying  deterministic  PDE  is 

du 


—  =  D\u  —  C.  \7  u  -h  aii,  t  >  0, 


where  u  =  Ut(x),  D  >  0,  V'  =  (V‘i _ V4)  and  a  are  constants 

'A  1 

dxi '  dx 


A  =  d  —  dimensional  Laplacian  and  V  —  f 


-V 


Let 


v; 


c,  =  p(x)  =  c  ^(•  >^1+ and  //„  =  L-(.V.  p(,r  )</x). 

The  cases  d  —  1,2  or  3  are  of  physical  interest. 

d  ~  1  (River  Pollution),  d  =  2  (or3)  (.Atmospheric  Pollution). 

We  impose  the  (Neumann)  boundary  conditions  (for  d  =  2), 

d  d  ()  d 

- — udO.x-i)  =  - — udf.x.))  =  - — udxj.O)  =  - — utixiJ)  =0  U  >  0). 
ox  I  ax  I  ox,  OXn 

The  problem  defines  a  positive,  self-adjoint  operator  denoted  by  —L  on  Ho  and  ( Tf  1 
is  the  semigroup  generated  by  L. 


Let  A;,  =  0.  A',  =  D(c"  -f  -mA)  i  — 


.d. 


~  (^)‘''‘e'-'^sin(  yiy -F  a 


J  ' 


w 


here  q)  =  tan  ‘(  — ;ff)-  J  >  I-  Then  A^, . and  Oj, . jj(x).  where 

'\i . 1  •  •  •  Oj|(.r,,). 

are  the  eigenvalues  and  eigenfunctions  of  —L.  {Oj,  j  J  is  a  CONS  in  Ho-  With  the 
help  of  these  we  <lefine  a  linear  space  ‘k  C  Ho  of  smootli  functions  in  Ho  and  we  have' 
a  chain 

<t>  C  . .  .  C  //p  C  .  .  .  C  Ho  C  //-I  C  . . .  C  C  .  ,  .  C  tl)'. 
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$  =  r\pHp\  H_p  =  Hp,  Hp  (Hilbert  space).  <I>  is  a  Frechel  space  which  is  also 
nuclear. 

The  natural  assumption  is  that  the  influx  of  pollution  follows  a  Poisson  process. 
More  precisely,  let  .V  be  a  Poisson  random  measure  on  R4.  x  x  R^.  (R.,.  =  (0.oc)) 
with  EN{dt  dx  da)  —  dt(j.{dx  da)  where  fx  is  the  intensity  measure  on  x  R+.  {At 
a  random  time  there  is  a  random  accretion  Aj  of  pollution). 

Then  the  stochastic  version  of  the  above  equation  is  given  by 

=  uo[v]  +  /  .4(s,  +  f  f  f  aip(x).\(ds  dx  da) 

Jo  Jo  Jx  JR+ 

(X  =  compensated  random  measure), 

.-1(5,  f’)[i^]  =  — +  /  /  aip{x)fx{dxda). 

Jx 

We  have  the  following  result  which  shows  that,  in  general,  the  SDE  has  a  solution  in 
so  that  we  do  not  have  to  seek  a  distributional  solution. 


Theorem  2  Suppose  E  ||  i/q  irwo<  and  /i  is  a  finite  measure.  Let 

1  /  /  a^p{x)fi{dx  da)\'  A-  [  /  a^Y{x)^(dx  da)  <  Const.  ||  ^  ||o  • 

Jx  JR+  Jx  Jr+ 

Then 

Ut  6  Ho  =  l^([0,f’]'^,pdx)  a.s.  Wt. 

For  d  =  1.  £  II  Uf  llo<  cc.  For  d  >  1.  £  ||  U(  ||o=  00  Vt  in  the  most  interesting  cases. 
Two  more  realistic  models  of  pollution  are  studied  (d=l); 


(i)  Pollution  emission  at  specific  '?ites  (Factories); 

(ii)  Pollution  model  with  an  upper  tolerance  level  for  chemical  concentration. 


(i)  Suppose  there  are  r  sites  A'l . 1\\  at  which  chemicals  are  deposited  in  terms 

of  independent  Poisson  streams  A',{t)  with  parameter  /,  >  0.  with  random  magnitudes 

=  1,2.. . .)  where  .4;  have  common  d.f.  F,(da)  arriving  in  the  vicinity  (A',  — 
e,.  A',  +  €,)  of  A',.  For  a  set  ,4  C  [0,/’]  and  B  C  R+  let 

r 

.V((0.i]x.4x£)  =  ^l.4(A'.)  £  ls(.40 

.=1  j=i 


with 

m(.4x£)  =  ^1,4(A'.)/.£.(£). 

1=1 

The  new  SDE  takes  the  form 


G{s.  ,  j,  ti  ).V((/.s  dx  da) 


where 


G(<,y,(x,a))[(^]  =  —  /  (l){y)dy  if  x  =  K„i  =  \,.. 

=  0  otherwise, 

-u[I^]  +  Qu[v?]  +  rn\ip], 

"  /.a, 


,r, 


A(i,y)[s?] 

m[sj] 


‘2e,  Jk,-u 


<p(y)dy 


and 


a,  =  /  aFi(da). 
'/R.+ 


(ii)  We  consider  a  simple  model  where  the  upper  level  for  pollution  is  a  fixed 
tolerance  function  ^(i).  Assume  that  the  change  of  chemical  concentration  does  not 
depend  on  the  locations  where  the  polluted  material  is  deposited.  We  then  have  a 
quasilinear  SDE 

Ut[v?]  =  +  [  {us[-L^]  +  auj[(^]}ds  +  /  /  a((fM  -  Us-[^])N{ds  da). 

Jo  Jo  JR+ 

N  being  a  Poisson  random  measure  with  intensity  measure  y,  on  R^.. 

If  the  initial  value  uq  is  smooth  and  a  <  0,  it  is  shown  that  the  above  quasi-linear 
SDE  has  a  solution  in  D([0,r],$)  and  furthermore,  under  suitable  conditions  that 
the  total  amount  of  pollution  cannot  exceed  a  prescribed  bound. 

Finally,  the  asymptotic  behavior  as  <  — »  oo  of  the  solution  of  the  SDE  in  (i)  is 
investigated  and  the  following  diffusion  approximation  is  obtained: 

.  . 

—  =  DAut  -  V  ■  s;^ut  +  Q!Ut+  Wtx 
at 

where  Wtx  is  space-time  Gaussian  white  noise 

where  Bj(t)  are  independent  real  Brownian  motions,  Fj(da)  and  is 

Dirac  measure  at  K. 


Research  on  interacting  Hilbert  space  valued  diffusions  carried  out  in  collaboration 
with  A.  Bhatt  and  R.L.  Karandikar  (Annual  Scientific  Report  1990-1991)  has  provided 
results  that  have  been  applied  to  the  asymptotic  behavior  of  interacting  neurons  in 
the  following  paper: 

4.  Stochastic  differential  equation  models  for  spatially  distributed  neurons 
and  propagation  of  chaos  for  interacting  systems  [4] 

Distribution  or  nuclear  space  valued  SDE’s  (diffusions  as  well  as  discontinuous 
equations)  are  discussed  as  stochastic  models  for  the  behavior  of  voltage  potentials 
of  spatially  distributed  neurons.  A  propagation  of  chaos  result  is  obtained  for  an 
interacting  system  of  Hilbert  space  valued  SDE’s. 


II.  APPLICATIONS  OF  NONLINEAR  WHITE  NOISE  STOCHASTIC 
ANALYSIS 

The  question  of  when  a  nonlinear  transformation  of  the  Wiener  measure  /i  is  ab¬ 
solutely  continuous  with  respect  to  ^  is  a  difficult  problem  that  has  been  outstanding 
since  the  time  of  Cameron  and  Martin  who  were  the  first  to  investigate  it.  The  most 
important  work  since  then  has  been  done  by  R.  Ramer  and  generalized  by  S.  Kusuoka. 

A  nonlinear  theory  of  white  noise  on  Hilbert  space  developed  by  R.L.  Karandikar 
and  myself  has  provided  a  new  way  to  approach  this  problem.  This  research  is 
presented  in  the  following  paper: 

5.  Nonlinear  transformations  of  the  canonical  Gauss  measure  on  Hilbert 
space  and  absolute  continuity  (with  R.L.  Karandikar)  [5] 

The  papers  of  R.  Ramer  (1974)  and  S.  Kusuoka  (1982)  investigate  conditions  under 
which  the  probability  measure  induced  by  a  nonlinear  transformation  on  abstract 
Wiener  space  (y,  H,  B)  is  absolutely  continuous  with  respect  to  the  abstract  Wiener 
measure  fi.  These  conditions  reveal  the  importance  of  the  underlying  Hilbert  space  H 
but  involve  the  space  B  in  an  essential  way.  The  present  paper  gives  conditions  solely 
based  on  H  and  takes  as  its  starting  point  a  nonlinear  transformation  T  =  I  +  F 
on  H.  New  sufficient  conditions  for  absolute  continuity  are  given  which  do  not  seem 
easily  comparable  with  those  of  Kusuoka  or  Ramer  but  are  more  general  than  those 
of  Buckdahn  (1991)  and  Enchev  (1991).  The  Ramer- Ito  integral  occurring  in  the 
expression  for  the  Radon- Nikodym  derivative  is  studied  in  some  detail  and,  in  the 
general  context  of  white  noise  theory,  it  is  shown  to  be  an  anticipative  stochastic 
integral  which,  under  a  stronger  condition  on  the  weak  Gateaux  derivative  of  F.  is 
directly  related  to  the  Ogawa  integral. 

III.  FEYNMAN  INTEGRALS;  FUNCTIONAL  INTEGRALS  OVER 
HILBERT  SPACES  RELATED  TO  THE  FEYNMAN  INTEGRAL. 


6.  Integration  over  Hilbert  spaces:  Examples  inspired  by  the  harmonic 
oscillator  (with  V.  Papanicolaou)  [6] 

The  work  is  joint  with  V.  Papanicolaou  and  is  briefly  described  under  his  heading. 

IV.  PREDICTION  THEORY  OF  SECOND  ORDER  STATIONARY 
RANDOM  FIELDS 

7.  Spectral  characterization  and  autoregressive  expansion  of  linear  pre¬ 
dictors  for  second  order  stationary  random  fields  (SOSRF),  Part  I  (with  J. 
Farshidi  and  V.  Mandrekar)  [7] 

The  problem  of  finding  spectral  criteria  for  autoregressive  (AR)  expansions  is  of 
great  practical  importance  for  single  parameter  stationary  time  series.  While  this 
problem  has  been  solved  satisfactorily  in  recent  years,  the  corresponding  problem  for 
SOSRF  has  only  now  come  to  the  forefront. 

The  major  difficulty  with  SOSRF  is  that  there  is  no  unique  definition  of  “past" 
and  “future”.  The  definitions  of  a  deterministic  and  purely  nondeterministic  random 
field  can  therefore  be  given  separately  for  the  horizontal,  vertical  and  “south  west" 
or  quarter  plane  past. 

An  AR  expansion  is  an  expansion  for  the  linear  least  squares  predictor  given  in 
terms  of  past  observations  rather  than  in  terms  of  an  innovation  sequence  based 
on  the  past.  The  original  .solution  of  the  Kolmogorov-Wiener  theory  is  based  on  the 


latter  and  is  less  easily  used  in  practiced  than  an  AR  expansion.  .A  paper  is  being 
prepared  on  the  work  done  to  date. 
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M.  ROSS  LEADBETTER 


Together  with  S.  Cambanis  and  G.  Kallianpur,  M.R.  Leadbetter  provided  contin¬ 
uing  direction  and  participation  in  the  research  activities  of  the  Center  for  Stochastic 
Processes.  Since  this  was  at  no  contract  cost  during  the  current  year,  a  brief  activity 
summary  is  given  here  rather  than  a  detailed  contract  reporting.  The  activities  are 
described  by  area  as  follows 


1.  Tail  inference  for  stochastic  sequences. 

Work  previously  reported  (cf.  Center  for  Stochastic  Processes  Technical  Report 
No.  292)  was  further  developed  jointly  with  H.  Rootzen.  This  concerns  problems  such 
as  estimation  of  parameters  of  exponentially  or  regularly  varying  tail  distributions, 
extremal  index,  tail  probabilities  and  quantiles.  This  work,  originally  planned  for  a 
single  publication,  is  now  being  expanded  into  two  parts. 

2.  Convergence  of  vector  random  measures. 

Research  with  S.  Nandagopalan  on  convergence  of  vector  random  measures  was  de¬ 
veloped  in  the  current  report  period  and  will  be  completed  in  the  subsequent  months. 
General  theorems  are  given  for  convergence,  with  particular  reference  to  the  random 
measure  formed  from  multilevel  exceedances  by  a  (nonstationary)  stochastic  process. 

3.  Processes  with  deterministic  peaks. 

Stationary  Gaussian  processes  have  the  property  that  high  peaks  have  an  in¬ 
creasingly  parabolic  asymptotic  form.  This  notion  can  be  generalized  leading  to  the 
concept  of  “deterministic  peaks”  -  where  the  time  above  a  high  level  (asymptotically) 
determines  that  above  any  higher  level.  This  ongoing  work  will  be  described  shortly 
in  a  paper  (joint  with  T.  Hsing). 


4.  Excursion  random  measures. 

Extreme  value  behavior  of  stochastic  sequences  can  be  summarized  by  limiting 
behavior  of  the  two  dimensional  point  process  formed  by  plotting  (a  suitably  normal¬ 
ized  version  of)  the  sequence  in  the  plane.  Substantial  effort  -  jointly  with  T.  Hsing  - 
has  been  put  into  the  development  of  a  corresponding  continuous  time  theory  in  this 
and  previous  reporting  periods  (cf.  Center  for  Stochastic  Processes  Technical  Report 
No.  350).  This  work  is  now  undergoing  revision  for  publication. 

5.  Applications. 

Work  (jointly  with  a  student,  L.S.  Huang)  was  initiated  on  the  application  of 
"exceedance  methods”  to  environmental  data.  Some  preliminary  time  series  modeling 
of  ozone  data  has  been  undertaken  and  will  be  continued  in  the  coming  year. 
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RAY  CHENG 


Professor  Ray  Cheng  of  the  Department  of  Mathematics  of  the  Univ'ersity  of 
Louisville  visited  the  Center  for  two  months  and  worked  on  the  structure  of  two- 
parameter  random  fields  which  is  relevant  to  the  problem  of  prediction.  He  completed 
the  following  technical  reports. 


1.  Outer  factorization  of  operator  valued  weight  functions  on  the  torus  [I] 

An  exact  criterion  is  derived  for  an  operator-valued  weight  function  W(e'\  e‘‘)  on 
the  torus  to  have  a  factorization  VV(e‘^,e“)  =  e“)'4>(e'*.  e'‘),  where  the  operator 

valued  Fourier  coefficients  of  <I>  vanish  outside  of  the  Helson-Lowdenslager  halfplane 
A  =  {(m,n)  G  ;  m  >  1}  U  {(0,n)  ;  n  >  0},  and  is  “outer"  in  a  related  sense. 
The  criterion  is  expressed  in  terms  of  a  regularity  condition  on  the  weighted  space 
L^{W)  of  vector  valued  functions  on  the  torus.  .A  logarithmic  integrability  test  is  also 
provided.  The  factor  $  is  explicitly  constructed  in  terms  of  Toeplitz  operators  and 
other  structures  associated  with  VV.  The  corresponding  version  of  Szegos  infimum  is 
given. 


2.  Operator  valued  functions  of  several  variables:  Factorization  and  in¬ 
variant  subspaces  [2] 

This  work  is  an  attempt  to  extend  the  classical  function  theory  on  the  Hardy  space 

to  certain  classes  of  operator  valued  functions  of  several  variables.  Of  course,  it  is 
impossible  to  carry  over  all  of  the  interesting  details.  Our  focus  is  to  adapt  the  notions 
of  inner  and  outer  functions,  so  as  to  preserve  two  basic  factorization  theorems.  We 
also  establish  a  sort  of  Beurling-Lax  theorem  to  describe  a  class  of  associated  invariant 
subspaces.  The  overall  approach  concerns  functions  on  the  torus,  which  generally 
cannot  be  realized  as  the  boundary  limits  of  analytic  functions  in  the  complex  sense. 
Accordingly,  our  techniques  are  chiefly  borrowed  from  multiple  Fourier  series  and 
shift  analysis. 
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ISSA  FAKHRE-ZAKERI 


Professor  Issa  Fakhre-Zakeri  of  the  Department  of  Mathematics  of  the  University 
of  Maryland  visited  the  cetner  during  the  1992  year.  He  worked  on  inference  problems 
for  stationary  linear  time  series  with  finite  variance  jointly  with  J.  Farshidi  [1.2)  and 
with  heavy  tails  jointly  with  S.  Cambanis  [3]. 

1.  A  central  limit  theorem  with  random  indices  for  stationary  linear  pro¬ 
cesses  [1] 

A  central  limit  theorem  with  random  indices  is  obtained  for  stationary  linear 
process  Xt  —  fJ.  =  where  {r)t)  are  independent  and  identically  distributed 

random  variables  with  mean  zero  and  finite  variance  and  \ 


2.  Limit  theorems  for  sample  covariances  of  stationary  linear  processes 
with  applications  to  sequential  estimation  [2] 

For  a  stationary  linear  process  the  strong  consistency  and  rate  of  convergence 
are  established  under  optimal  conditions  for  the  asymptotic  variance  of  their  sample 
mean.  Applications  are  made  to  the  problem  of  sequential  point  and  fixed  width 
confidence  interval  estimation  of  the  mean  of  a  stationary  linear  process. 

3.  On  prediction  of  heavy-tailed  autoregressive  sequences:  Regression 
versus  best  linear  prediction  [3] 

The  prediction  of  heavy-tailed  first  order  autoregressive  sequences  is  considered. 
In  forward  time  the  regression  on  all  p«ist  values  is  the  same  as  the  one-step  regression 
on  the  previous  value,  which  is  in  fact  linear.  In  reversed  time  the  regression  on  all 
future  values  is  the  same  as  the  one-step  regression  on  the  immediate  future  value 
(i.e.  the  Markovian  property  is  retained)  and  we  show  that  it  is  linear  if  and  only 
if  the  innovations  have  a  semistable  distribution.  This  answers  a  question  posed  by 
Rosenblatt  (1992)  who  considered  sequences  with  finite  second  moment  and  showed 
that  regression  with  time  reversed  is  linear  if  and  only  if  the  innovations  are  Gaussian. 

When  the  distribution  of  the  innovations  is  non-Gaussian  stable,  then  both  re¬ 
gressions  in  forward  and  reversed  time  are  linear,  but  while  the  forward  regression  is 
the  best  linear  predictor,  the  regression  with  time  reversed  is  not!  The  performance 
of  linear  regression  predictor  is  compared  in  this  case  with  that  of  the  best  linear 
predictor. 
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JAMSHID  FARSHIDI 


Dr.  Jamshid  Farshidi  from  the  Department  of  Probability  and  Statistics  of  Michi¬ 
gan  State  Unviersity  spent  the  academic  year  cis  a  postdoctoral  visitor  to  the  center. 
He  worked  on  the  problem  of  prediction  of  stationary  time  series  [1]  and  of  random 
fields  [2]  and  on  inference  for  stationary  linear  time  series  (jointly  with  I.  Fakhre- 
Zakeri)  [3,4].  He  has  also  begun  working  on  heavy  tailed  stationary  time  series  and 
more  specifically  on  the  prediction  of  harmonizable  stationary  stable  processes. 

1.  Autoregressive  expansion  of  the  linear  predictor  for  stationary  stochas¬ 
tic  processes  [1] 

The  principal  problems  considered  are  the  existence  and  uniqueness  of  an  au¬ 
toregressive  expansion  of  the  linear  predictor  for  a  discrete  stationary  process  with 
spectral  density  /  and  optimal  factor  and  the  invertibility  of  the  process  .V.  The 
main  results  are: 

(1)  the  equivalence  of  the  strong  convergence  of  an  autoregresive  series  to  the 
linear  predictor,  with  its  boundedness,  and  with  its  weak  convergence; 

(2)  the  uniqueness  of  an  autoregressive  expansion; 

(3)  the  equivalence  of  an  autoregressive  expansion  with  the  invertibility  of  the 
process: 

(4)  the  sufficiency  of  the  condition  (1//)  G  for  the  existence,  convergence, 
uniqueness  of  the  autoregressive  expansion  and  the  invertibility  of  the  process; 

(5)  a  necessary  condition  based  on  and  /  for  the  existence,  uniqueness,  and 
convergence  of  an  autoregressive  expansion,  and  invertibility  of  the  process. 

2.  Spectral  characterization  and  autoregressive  expansion  of  linear  pre¬ 
dictors  for  second  order  stationary  random  fields  (SOSRF),  Part  I  (with  G. 
Kallianpur  and  V.  Mandrekar)  [2] 

The  problem  of  finding  spectral  criteria  for  autoregressive  (.AR)  expansions  is  of 
great  practical  importance  for  single  parameter  stationary  time  series.  While  this 
problem  has  been  solved  satisfactorily  in  recent  years,  the  corresponding  problem  for 
SOSRF  has  only  now  come  to  the  forefront. 

The  major  difficulty  with  SOSRF  is  that  there  is  no  unique  definition  of  "past" 
and  "future".  The  definitions  of  a  deterministic  and  purely  nondeterministic  random 
field  can  therefore  be  given  separately  for  the  horizontal,  vertical  and  "south  west" 
or  quarter  plane  past. 

An  AR  e.xpansion  is  an  expansion  for  the  linear  least  squares  predictor  given  in 
terms  of  past  observations  rather  than  in  terms  of  an  innovation  sequence  based 
on  the  past.  The  original  solution  of  the  Kolmogorov-Wiener  theory  is  based  on  the 
latter  and  is  less  easily  used  in  practiced  than  an  .AR  expansion.  .\  paper  is  being 
prepared  on  the  work  done  to  date. 

3.  A  central  limit  theorem  with  random  indices  for  stationary  linear  pro¬ 
cesses  [3] 

A  central  limit  theorem  with  random  indices  is  obtained  for  stationary  linear 
process  Xt  —  where  {r/(}  are  independent  and  identically  distributed 

random  variables  with  mean  zero  and  finite  variance  and  |nj|  <  dc. 

4.  Limit  theorems  for  sample  covariances  of  stationary  linear  processes 
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with  applications  to  sequential  estimation  [4] 

The  strong  consistency  and  rate  of  convergence  are  established  under  optimal 
conditions  for  the  asymptotic  variance  of  the  sample  mean  of  a  stationary  linear 
process  .  Applications  are  made  to  the  problem  of  sequential  point  and  fixed  width 
confidence  interval  estimation  of  the  mean  of  a  stationary  linear  process. 
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LARS  HOLST 


Professor  Lars  Holst  of  the  Royal  Institute  of  Technology,  in  Stockholm,  visited 
the  Center  in  January  1992.  Professor  Holst  is  an  authority  on  Poisson  approxi¬ 
mations  and  interacted  with  faculty  and  visitors,  also  presenting  a  colloquium  on 
approximation  of  Stein’s  Method. 
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CHRISTIAN  HOUDRE 


Dr.  Christian  Houdre  of  the  Department  of  Mathematics  of  the  University  of 
Maryland,  now  at  the  Department  of  Statistics  of  Stanford  University,  visited  the 
Center  for  two  months.  He  worked  primarily  on  the  ramification  of  wavelets  in 
stochastic  processes  [1,2,3],  the  latter  being  ongoing  collaboration  with  S.  Camba- 
nis.  He  also  worked  on  stable  stochastic  processes  jointly  with  M.  Hernandez  [4]  and 
on  variance  inequalities  for  functions  of  Gaussian  random  variables  jointly  with  A. 
Kagan  [5]. 

1.  Wavelets,  probability  and  statistics:  Some  bridges  [1] 

The  role  of  some  wavelet  methods  in  probability  and  statistics  is  illustrated  via 
a  sample  of  three  problems:  We  show  how  properties  of  processes  can  be  read  off 
properties  of  their  wavelet  transform.  We  discuss  how  the  missing  data  problem  can 
be  approached  via  frames  of  complex  exponentials.  We  explain  how  wavelets  can  be 
used  to  span  classes  of  admissible  estimators  in  non-parametric  function  estimation. 
It  is  also  the  purpose  of  this  paper  to  show  that  bridges  can  be  crossed  in  the  other 
direction.  Random  products  of  matrices  determine  the  smoothness  of  compactly 
supported  wavelets.  Non  stationary  prediction  theory  gives  new  results  on  frames  in 
Hilbert  space. 

2.  Path  reconstruction  of  processes  from  missing  and  irregular  samples  [2] 

A  criterion  is  provided  for  the  reconstruction  of  the  paths  of  non-stationary  band- 
limited  processes  using  irregularly  spaced  samples  by  means  of  an  interpolation  for¬ 
mula.  Its  rate  of  convergence  is  studied  along  with  its  truncation  error.  These  results 
provide  irregular  sampling  theorems  for,  say,  deterministic  signals  corrupted  by  addi¬ 
tive  noise,  and  a  potential  solution  to  the  missing  data  problem:  interpolation  from 
sparse  or  missing  data  can  be  achieved  under  a  density  condition.  The  analysis  in¬ 
volves  classical  results  on  non-harmonic  Fourier  series  as  well  as  more  recent  results 
on  frames  and  wavelets. 


3.  Wavelet  transforms  of  random  processes  [3] 

A  study  has  been  initiated  of  the  properties  of  wavelet  transforms  of  random 
processes  whose  sampled  values  appear  as  coefficients  in  the  wavelet  approximation 
of  the  process  at  a  given  resolution.  .A  natural  question  is  which  properties  of  the 
process  are  inherited  to  its  wavelet  transform,  and,  conversely,  which  properties  of  the 
process  can  be  read-off  properties  of  its  wavelet  transform.  For  random  processes  with 
finite  second  moment,  properties  such  as  periodicity,  stationarity.  harmonizability. 
and  self-similarity,  are  characterized  by  means  of  analogous  properties  of  their  wavelet 
transforms  at  some  scale:  The  properties  of  the  wavelet  transform  characterize  the 
corresponding  properties  of  the  increments  of  the  process  of  order  equal  to  the  order 
of  regularity  of  the  analyzing  wavelet. 

4.  Disjointness  results  for  some  classes  of  stable  processes  [4] 

The  disjointness  of  two  classes  of  stable  stochastic  processes:  moving  averages  and 
Fourier  transforms  is  discussed.  Results  on  the  incompatibility  of  these  two  repre¬ 
sentations  date  back  to  Urbanik  (1964).  Here  we  extend  varous  earlier  disjointness 
results  to  encompass  larger  classes  of  processes,  allowing  e.g.  the  noise  of  a  moving 
average  process  to  be  nonstationary  and  showing  that  all  moving  average  processes 
are  Fourier  transforms  in  the  summability  sense. 
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5.  Variance  inequalities  for  functions  of  Gaussian  variables  [5] 

When  X  is  a  standard  Gaussian  random  variable  and  G  an  absolutely  continuous 
function,  the  inequality  Var[G(X)]  <  E[G'{X)Y  was  proved  in  Nash  (1958)  and  later 
rediscovered  in  Brascamp  and  Lieb  (1976)  as  a  special  case  of  a  general  inequality 
in  ChernofF  (1981).  All  the  proofs  are  based  on  properties  of  the  Gaussian  density. 
By  using  the  characteristic  function  rather  than  the  density,  generalizations  with 
higher  orcier  derivatives  are  obtained.  The  method  also  establishes  potentially  useful 
connections  with  Karlin's  total  positivity. 
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TAILEN  HSING 


Professor  Tailen  Hsing  of  the  Statistics  Department  of  Texas  A  &  M  University 
visited  the  Center  for  the  1990-91  academic  year.  In  addition  to  his  work  reported 
in  last  year’s  annual  report,  in  the  following  paper  he  extended  and  completed  the 
work  on  the  estimation  of  the  spectral  density  of  harmonizable  stable  processes  of 
Cambanis  and  Masry  (1984). 


1.  Limit  theorems  for  stable  processes  with  application  to  spectral  density 
estimation  [1] 

It  is  shown  that  for  a  nearly  stationary  moving  average  a-stable  process  V’  and  for 
each  fixed  0  <  p  <  oo,  a  weighted  average  of  |K(f)j'’ over  [—7’,  T]  has  an  asymptotically 
(2  A  ^)-stable  distribution  as  T  — »  oo.  This  is  a  partial  extension  of  the  limit  theorems 

considered  in  Davis  (1983)  and  LePage,  Woodroofe  and  Zinn  (1981).  Applications  of 
the  results  are  made  in  the  context  of  spectral  density  estimation  of  a  harmonizable 
a-stable  process.  The  spectral  density  estimator  is  the  smoothed  v'crsion  of  the  pth 
absolute  power  of  the  tapered  Fourier  transform  proposed  in  Cambanis  and  Masry 
(1984)  and  proven  consistent  when  0  <  p  <  a/2.  Here  its  asymptotic  distribution 
is  derived  and  is  shown  to  be  normal  when  0  <  p  <  a/2  and  (a/p)-stable  when 
a/2  <  p  <  a.  Also  the  best  possible  rates  of  convergence  are  determined  and  show 
that  the  rate  of  convergence  is  faster  for  p  in  (0,a/2). 
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HARRY  HURD  -  ANDRZEJ  RUSSEK 


Dr.  Harry  Hurd  continued  the  systematic  study  of  non-stationary  processes  which 
are  periodically  correlated  jointly  with  Dr.  Andrzej  Russek  of  the  Polish  Academy  of 
Science,  Sopot,  who  has  been  a  visitor  to  the  Center  since  January  1992. 


1,  Stepanov  almost  periodically  correlated  and  almost  periodically  unitary 
processes  [1] 

VVe  extend  the  structure  and  properties  of  almost  periodically  correlated  (APC) 
and  almost  periodically  unitary  (APU)  processes,  which  were  defined  in  the  sense  of 
Bohr,  to  a  larger  class  of  processes  for  which  the  sense  of  almost  periodicity  is  that 
of  Stepanov.  These  processes  are  not  necessarily  continuous  in  quadratic  mean,  as 
are  the  Bohr  APC  and  APU  processes,  but  yet  exhibit  a  sense  of  almost  periodic¬ 
ity.  For  example,  processes  formed  by  amplitude  modulation  f{t)X{t)  or  time-scale 
modulation  X{t  /(<))  of  a  wide  sense  stationary  process  X(t)  by  a  Stepanov  .APU 
and  APC.  The  principal  results  on  APC  and  APU  processes  are  extended  to  the 
new  class.  We  extend  Gladyshev’s  characterization  of  APC  correlation  functions  to 
Stepanov  APC  processes  and  show  that  their  correlation  functions  are  completely  rep¬ 
resented  by  a  Fourier  series  having  a  countable  number  of  coefficient  functions  that 
are  Fourier  transforms  of  complex  measures.  We  show  that  Stepanov  .APU  processes 
are  also  Stepanov  APC  and  are  given  by  X(t)  =  U{t)[P(t)]  where  {U{t).  <  6  R}  is  a 
strongly  continuous  group  of  unitary  operators  and  P(t)  is  a  vector- valued  Stepanov 
almost  periodic  function.  As  in  the  case  of  Bohr  APU  processes,  the  preceeding  fact 
leads  to  representations  of  X(t)  based  on  the  spectral  theory  for  unitary  operators 
and  for  Stepanov  almost  periodic  functions. 


2.  Almost  periodically  correlated  processes  on  LCA  groups  [2] 

For  an  almost  periodic  covariance  R{t  +  r.t)  =  E{X{t  T),V(t)}  of  a  second 
order  stochastic  process  X{t)  indexed  by  an  LC.A  group  G.  we  show  that  the  means 

a(/\,r)  =  +  r.  <)/\(t)}  are  Fourier  transforms  of  signed  measures  with  finite 

total  variation.  We  examine  conditions  under  which  X{t)  or.  more  precisely,  its 
correlation  R{t  +  T,t),  has  a  countable  set  of  spectral  characteristic  exponents  (or  fre¬ 
quencies).  We  also  consider  the  problem  of  finding  a  stationarizing  shift  and  exhibit 
a  class  of  G-valued  random  variables  9  such  that  Y'{t)  =  X(t  +  9)  is  stationary.  Fi¬ 
nally  we  characterize  the  almost  periodically  correlated  processes  among  the  stronglv 
harmonizable  ones. 
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RAJEEVA  L.  KARANDIKAR 


Professor  Karandikar  of  the  Indian  Statistical  Institute,  Delhi,  visited  the  Center 
for  four  months  in  1992.  In  addition  to  the  completion  of  work  done  jointly  with  A. 
Bhatt  and  G.  Kallianpur  [1]  and  partially  described  in  the  last  report,  he  collaborated 
with  V.G.  Kulkarni  of  the  Operations  Research  department  on  the  study  of  a  second- 
order  fluid  flow  model  [2]  and  with  G.  Kallianpur  on  nonlinear  transformations  of 
abstract  Wiener  measure  [3]. 


1.  On  interacting  systems  of  Hilbert  space  valued  diffusions  [1] 

.A  nonlinear  Hilbert  space  valued  stochcistic  differential  equation  where  [L 
being  the  generator  of  the  evolution  semigroup)  is  not  nuclear  is  investigated  in  the 
paper.  Under  the  assumption  of  nuclearity  of  the  existence  of  a  unique  solution 
lying  in  the  Hilbert  space  H  has  been  shown  by  Dawson  in  an  early  paper.  When  L~^ 
is  not  nuclear,  a  solution  in  most  cases  lies  not  in  H  but  in  a  larger  Hilbert,  Banach  or 
nuclear  space.  Part  of  the  motivation  of  the  present  paper  is  to  prove  under  suitable 
conditions  that  a  unique  strong  solution  can  still  be  found  to  lie  in  the  space  H  itself. 
Uniqueness  of  the  weak  solution  is  proved  without  moment  assumptions  on  the  initial 
random  variable. 

.A  second  problem  considered  is  the  asymptotic  behavior  of  the  sequence  of  em¬ 
pirical  measures  determined  by  the  solutions  of  an  interacting  system  of  //-valued 
diffusions.  It  is  shown  that  the  sequence  converges  in  probability  to  the  unique 
so'ution  .\o  of  the  martingale  problem  posed  by  the  corresponding  McKean- Vlasov 
equation. 


2.  Second-order  fluid  flow  model  of  a  data-buffer  ir*  random  environment 

[21 

This  paper  considers  a  stochastic  model  ot  a  data-buffer  in  a  telecommunication 
network.  Let  X{t)  be  the  buffer-content  at  time  t.  The  {X(t).  t  >  0}  process 
depends  on  a  finite  state  continuous  time  Markov  process  ^  >  0}  as  follows: 

during  the  time-intervals  when  Z{t)  is  in  state  i,  X{t)  is  a  Brownian  motion  with 
drift  /i,,  variance  parameter  erf  and  a  reflecting  boundary  at  zero.  This  paper  studies 
the  steady  state  analysis  of  the  bivariate  process  {(,Y(f).  Z{t)),  i  >  0}  in  terms  of  the 
eigenvalues  and  eigenvectors  of  a  non-linear  matrix  system.  Algorithms  are  developed 
to  compute  the  steady  state  distributions  as  well  as  moments. 

Numerical  work  is  reported  to  show  that  the  variance  parameter  has  a  dramatic 
effect  on  the  buffer  content  process.  Thus  buffer  sizing  done  with  first  order  fluid  flow 
models  (with  zero  variance  parameters)  should  be  used  with  care. 

3.  Nonlinear  transformations  of  the  canonical  Gauss  measure  on  Hilbert 
space  and  absolute  continuity  [3j 

The  papers  of  R.  Ramer  ( 1974)  and  S.  Kusuoka  ( 1982)  investigate  conditions  under 
which  the  probability  measure  induced  by  a  nonlinear  transformation  on  abstract 
Wiener  space  (7,  //,  B)  is  absolutely  continuous  with  respect  to  the  abstract  Wiener 
measure  /i.  These  conditions  reveal  the  importance  of  the  underlying  Hilbert  space  H 
but  involve  the  space  B  in  an  essential  way.  The  present  paper  gives  conditions  solely 
based  on  H  and  takes  as  its  starting  point  a  nonlinear  transformation  T  —  I  +  F 
on  H .  New  sufficient  conditions  for  absolute  continuity  are  given  which  do  not  seem 
easily  comparable  with  those  of  Kusuoka  or  Ramer  but  are  more  general  than  those 
of  Buckdahn  (1991)  and  Enchev  (1991).  The  Ramer-Ito  integral  occurring  in  the 
expression  for  the  Radon-Nikodym  derivative  is  studied  in  some  detail  and.  in  the 
general  context  of  white  noise  theory,  it  is  shown  to  be  an  anticipative  stochastic 
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integral  which,  under  a  stronger  condition  on  the  weak  Gateaux  derivative  of  F,  is 
directly  related  to  the  Ogawa  integral. 
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JEAN-CLAUDE  MASSE 


Professor  Jean-Claude  Masse  of  the  Department  of  Mathematics  and  Statistics  of 
the  Unviersity  of  Laval  visited  the  Center  for  six  weeks  and  studied,  jointly  with  C.A. 
Leon,  the  properties  of  the  simplicial  median  of  Oja  with  a  view  towards  the  study 
of  medians  of  random  processes. 

1.  La  mediane  simpliciale  d’Oja:  existence,  unicite  et  stabilite  [1] 

Oja  (1983)  examined  various  ways  of  measuring  location,  scatter,  skewness  and 
kurtosis  for  multivariate  distributions.  Among  other  measures  of  location,  he  intro¬ 
duced  a  generalised  median,  the  Oja  median.  We  study  three  fundamental  theoretical 
properties  of  that  median:  existence,  uniqueness  and  consistency. 
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DITLEV  MONRAD 


Professor  Ditlev  Monrad  of  the  Department  of  Statistics  of  the  University  of  Illi¬ 
nois  visited  the  Center  for  four  months  and  worked  on  sample  function  properties 
of  Levy  processes  with  local  time,  and  of  fractional  Brownian  motion  and  other  lo¬ 
cally  nondeterministic  Gaussian  processes.  The  latter  work  was  done  jointly  with  H. 
Rootzen. 


1.  Some  uniform  dimension  theorems  for  the  sample  functions  of  Levy 
processes  with  local  times  [1] 

When  A  is  a  real-valued,  strictly  stable  Levy  process  of  index  a,  1  <  a  <  2,  we 
show  that  with  probability  one, 


dim  A-'(F)  =  1  -  -  -f  -  dim  (F), 
a  Of 

simultaneously  for  all  real  Borel  sets  F.  The  result  is  also  extended  to  general  real¬ 
valued  Levy  processes  with  local  time. 

2.  Small  values  of  fractional  Brownian  motion  and  locally  nondeterministic 
Gaussian  processes  [2] 

A  centered  Gaussian  process  {5^(0  •  ^  ^  0}  with  covariance  function  proportional 
to  |{|3|“  +  1^1"  —  l-s  —  for  a  €  (0, 2)  is  called  a  fractional  Brownian  motion  (fBm). 
Chung  type  laws  of  the  iterated  logarithm  are  proved  for  fBm’s:  for  any  sample  path, 
there  are  arbitrarily  large  values  of  t  for  which  (Ba(s)  :  0  <  s  <  t}  is  confined  to  the 
interval  ±  const.  <"^^(log  log  but  this  is  not  true  for  any  narrower  intervals. 

A  corresponding  result  holds  for  small  values  of  L  Let  A/(C  =  maxo<s<tlBa(s)l.  For 
the  proof  the  following  bounds,  valid  for  small  c’s,  and  constants  0  <  c  <  C  are  found 

^ 

for  the  probability  that  the  process  is  flat.  They  hold  for  strongly  locally  nondeter¬ 
ministic  Gaussian  processes  whose  incremental  variances  over  intervals  of  length  h 
are  roughly  proportional  to  /i". 
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JOHN  P.  NOLAN 


Professor  John  Nolan  of  the  Department  of  Mathematics  and  Statistics  of  the 
American  University  in  Washington,  DC,  completed  some  work  on  multidimensional 
stable  distributons  which  Wcis  substantially  performed  during  a  1990  visit  to  the 
Center  but  was  not  included  in  the  1989-90  annual  report. 

1.  Approximation  of  multidimensional  stable  densities  [1] 

Stable  densities  in  two  or  more  variables  do  not  generally  have  explicit  formula. 
One  way  of  characterizing  these  distributions  is  by  a  spectral  measure.  Our  main 
result  shows  that  densities  and  probabilities  can  be  uniformly  approximated  by  ap¬ 
proximating  the  spectral  measure  with  a  discrete  spectral  measure  having  a  finite 
number  of  atoms.  A  concrete  formula  is  given  for  the  number  of  atoms  needed  and 
their  weights,  this  can  be  used  to  numerically  calculate  multidimensional  stable  den¬ 
sities.  Sample  graphs  of  two  dimensional  stable  densities  with  dependence  are  given. 
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TOMMY  NORBERG 


Dr.  Tommy  Norberg  of  the  University  of  Goteborg  visited  for  a  one  month  period 
in  November  1991,  primarily  to  collaborate  with  M.R.  Leadbetter  on  problems  in 
point  process  theory.  Dr.  Norberg  is  a  foremost  authority  on  the  theory  of  random 
sets,  which  provide  a  useful  alternative  framework  complementing  those  of  point 
processes  and  random  measures. 

During  the  visit  Dr.  Norberg  and  M.R.  Leadbetter  worked  together  in  three  areas; 

(a)  Foundations  for  a  non  topological  (or  minimally  topological)  theory  of  point  pro¬ 
cesses 

(b)  The  potential  use  of  random  sets  in  applied  areas  -  such  cis  minefield  modeling 
in  defense  applications 

(c)  Planning  for  a  volume  on  point  processes,  random  sets  and  random  measures.  This 
will  describe  the  different  structural  frameworks,  their  relationships,  the  usefulness  of 
each  view  and  some  of  their  applications. 


Work  has  continued  in  these  areas  since  the  visit  and  a  start  has  been  made  on 
the  writing  under  (c). 


JONNY  OLSSON 


Dr.  Jonny  Olsson  (University  of  Lund,  Sweden)  weis  a  junior  visitor,  supported 
from  Swedish  sources  for  November  1991.  He  collaborated  with  H.  Rootzen  on  Markov- 
random  fields  for  vision. 


1.  Image  Modeling  [1] 

Markov  random  field  models  for  vision  were  developed  jointly  with  J.  Olsson, 
a  junior  short  term  visitor  (noted  above),  describing  the  theory  and  application  to 
peripheral  vision  assessment,  and  is  summarized  as  follows. 

Measurement  of  the  patient’s  “seeing  threshold”  at  different  points  in  the  visual 
field  is  an  important  diagnostic  tool  for  glaucoma  and  other  diseases.  A  Markov 
random  field  model  is  developed  and  used  for  efficient  estimation  of  the  thresholds 
and  simultaneously  for  classification  of  the  measured  points  as  “normal”  or  “defec¬ 
tive”.  The  model  allows  for  nonhomogeneous  spatial  dependence  and  nonsymmetric 
marginal  distributions  and  hcis  physically  interpretable  parameters.  Maximum  a  pos¬ 
teriori  threshold  estimation  of  visual  fields  results  in  13%  -  31%  reduction  of  mean 
square  error  (depending  on  the  patient  population)  as  compared  to  currently  used 
procedures  and  in  a  fair  agreement  between  true  and  estimated  defect  status. 

Nonstandard  features  of  the  problem  are:  (i)  the  picture  is  small,  (ii)  there  is 
a  nonhomogeneous  directional  dependence,  and  (iii)  thresholds  are  only  measured 
indirectly,  by  binary  responses  to  questions,  where  the  probability  of  response  depends 
on  the  threshold  and  the  stimulus  level. 
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VASSILIS  PAPANICOLAOU 


Professor  Papanicolaou  of  the  Department  of  Mathematics  of  Duke  University 
[and  now  at  the  Department  of  Mathematics  of  Wichita  State  University]  was  at  the 
Center  for  two  months  in  the  summer  of  1992.  His  interest  in  Feynman  integrals 
led  to  the  study  of  some  problems  of  integration  over  Hilbert  space  and  extensions 
involving  multiparameter  Gaussian  processes  of  previous  work  by  G.  Kallianpur.  D. 
Kannan  and  R.L.  Karandikar  (Analytic  and  sequential  Feynman  integrals  on  abstract 
Wiener  and  Hilbert  spaces,  and  a  Cameron-Martin  formula,  Ann.  Inst.  H.  Poincare, 
21,  1985,  323-361). 

1.  Integration  over  Hilbert  spaces:  Examples  inspired  by  the  harmonic 
oscillator  [1] 

The  research  produced  some  examples  of  functional  integrals  over  Hilbert  spaces 
where  the  integrand  is  analogous  to  the  one  for  the  quantum  mechanical  harmonic  os¬ 
cillator.  In  one  czise  the  continuum  limit  of  a  sequence  of  coupled  harmonic  oscillators 
is  considered. 
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ROLAND  PERFEKT 


Dr.  Rx)land  Perfekt  (University  of  Lund,  Sweden)  was  a  junior  visitor,  supported 
from  Swedish  sources  for  the  period  September-October  1991.  He  collaborated  with 

H.  Rootzen  on  extremal  properties  of  stationary  Markov  chains. 

I,  Extremal  behaviour  of  stationary  Markov  chains  with  applications  [1] 

Extremal  behaviour  of  real-valued,  stationary  Markov  chains  is  studied  under 
rather  general  assumptions.  Conditions  are  obtained  for  convergence  in  distribution  of 
multi-level  exceedance  point  processes  associated  with  suitable  families  of  ‘increasing 
levels’.  Although  applicable  to  general  stationary  sequences,  these  conditions  are 
tailored  for  Markov  chains  and  are  seen  to  hold  for  a  large  class  of  chains.  The 
extra  assumptions  needed  are  that  the  marginal  distributions  belong  to  the  domain 
of  attraction  of  some  extreme  value  law  together  with  rather  weak  conditions  on  the 
transition  probabilities.  Also,  a  complete  convergence  result  is  given.  The  results 
are  applied  to  a  discrete-time  Lindley  process,  to  an  A/?(l)  process  with  uniform 
margins  and  to  solutions  of  a  first  order  stochastic  difference  equation  with  random 
coefficients. 
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HOLGER  ROOTZEN 


Professor  Holger  Rootzen  (University  of  Lund,  Sweden)  spent  12  months  as  a 
senior  visitor  to  the  Center.  In  this  period  he  also  arranged  shorter  visits  by  two  junior 
Swedish  statisticians  J.  Olsson  and  R.  Perfekt  (supported  from  Swedish  sources). 

Professor  Rootzen's  activities  were  divided  into  five  areas  as  follows. 


1.  Tail  estimation  for  stationary  sequences  [1] 

Holger  Rootzen  and  M.R.  Leadbetter  collaborated  in  extending  work  reported 
previously  (as  CSP  Report  292  with  L.  de  Haan)  on  the  estimation  of  parameters 
associated  with  high  values  of  stochastic  sequences.  These  include  the  “extremal 
index'’,  the  parameter  of  a  regularly  varying  tail  distribution,  tail  probabilities  and 
quantiles,  under  dependence  conditions  such  as  strong  mixing. 

It  is  planned  that  this  work  will  be  reported  in  two  papers  (currently  under  prepa¬ 
ration),  dividing  and  extending  that  in  the  CSP  Report  292. 

2.  Extremal  properties  of  Markov  chains  [2] 

Professor  Rootzen  worked  with  R.  Perfekt  on  extremal  properties  of  Markov 
chains.  This  is  described  in  the  research  activity  summary  for  R.  Perfekt. 


3.  Fractional  Brownian  motion  [3]  (with  D.  Monrad) 

Joint  work  was  conducted  on  sample  function  properties  of  fractional  Brownian 
motion  and  locally  nondeterministic  Gaussian  processes.  This  work  is  described  in 
CSP  Tech  Report  No.  361.  whose  contents  are  summarized  as  follows: 

A  centered  Gaussian  process  {Ba{t)  :  <  >  0}  with  covariance  function  proportional 
to  5{|sl“  -f  It]®  —  |3  —  tl“}  for  a  €  (0,2)  is  called  a  fractional  Brownian  motion  (fBm). 
Chung  type  laws  of  the  iterated  logarithm  are  proved  for  fBm's:  for  any  sample  path, 
there  are  arbitrarily  large  values  of  t  for  which  {Ba{s)  :  0  <  s  <  t)  is  confined  to  the 
interval  ±  const.  f‘*/^(Iog  log  but  this  is  not  true  for  any  narrower  intervals. 

A  corresponding  result  holds  for  small  values  of  t.  Let  ,V/(<)  =  maxo<,<tl5a(s)|.  For 
the  proof  the  following  bounds,  valid  for  .small  c’s,  and  constants  0  <  c  ^  G  are  found 

g-cu-^/o  <  p(M(t)  <t)< 

for  the  probability  that  the  process  is  flat.  They  hold  for  strongly  locally  nondeter¬ 
ministic  Gaussian  processes  whose  incremental  variances  over  intervals  of  length  h 
are  roughly  proportional  to  h'^ . 

4.  Image  Modeling  [4] 

Markov  random  field  models  for  vision  were  developed  jointly  with  J.  Olsson. 
a  junior  short  term  visitor  (noted  above),  describing  the  theory  and  application  to 
peripheral  vision  assessment,  and  is  summarized  as  follows. 

Measurement  of  the  patient’s  “seeing  threshold”  at  different  points  in  the  visual 
field  is  an  important  diagnostic  tool  for  glaucoma  and  other  diseases.  .A  Markov 
random  field  model  is  developed  and  used  for  efficient  estimation  of  the  thresholds 
and  simultaneously  for  classification  of  the  meeisured  points  as  "normal”  or  "defec¬ 
tive”.  The  model  allows  for  nonhomogeneous  spatial  dependence  and  nonsymmetric 
marginal  distributions  and  has  physically  interpretable  parameters.  Maximum  a  pos¬ 
teriori  threshold  estimation  of  visual  fields  results  in  13%  -  31%  reduction  of  mean 
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square  error  (depending  on  the  patient  population)  as  compared  to  currently  used 
procedures  and  in  a  fair  agreement  between  true  and  estimated  defect  status. 

Nonstandard  features  of  the  problem  are:  (i)  the  picture  is  small,  (ii)  there  is 
a  nonhomogeneous  directional  dependence,  and  (iii)  thresholds  are  only  measured 
indirectly,  by  binary  responses  to  questions,  where  the  probability  of  response  depends 
on  the  threshold  and  the  stimulus  level. 


5.  Related  statistical  questions 

Professor  Rootzen  conducted  related  statistical  research  in 

(i)  Quantile  estimation  in  a  nonparametric  component  of  variance  frame¬ 
work  with  applications  to  vision  problems  [5]. 

(ii)  Proportional  hazard  testing  related  to  strength  of  materials  [6]. 

His  results  under  (i)  are  summarized  as  follows: 

Quantile  estimators  for  a  non- parametric  components  of  variance  situation  are 
proposed  and  consistency  and  asymptotic  normality  is  proved.  Situations  with  differ¬ 
ent  numbers  of  measurements  for  different  subjects  are  considered.  Measurements  on 
separate  subjects  are  assumed  independent  while  measurements  on  the  same  subject 
have  a  fixed  dependence.  The  estimators  are  obtained  by  inverting  weighted  empiri¬ 
cal  distribution  functions.  An  "optimal”  estimator  is  derived  by  choosing  weights  to 
minimize  the  variance  of  the  weighted  empirical  distribution  function.  The  resulting 
weights  depend  on  unknown  parameters.  However,  these  weights  may  be  estimated 
from  data  without  affecting  asymptotic  performance.  A  simple  estimator  based  on 
within  subject  averages  i  investigated.  Small  sample  properties  are  studied  by 
simulation,  and  as  an  illustration  the  estimators  are  applied  to  normal  limits  for 
differential  light  sensitivity  of  the  eye. 

The  work  on  proportional  hazards  (joint  with  A.  Deis)  provided  a  A'-sample  test 
for  proportional  hazards  and  is  described  in  detail  as  follows: 

A  test  (or  proportionality  of  the  cumulative  hazard  functions  in  k  >  2,  possibly 
censored,  samples  is  proposed.  The  test  does  not  use  dummy  time-dependent  co¬ 
variates  or  partitions  of  the  time  axis.  It  extends  a  test  of  Wei  (1984)  from  2  to  k 
samples,  and  for  k  =  2  gives  an  alternative  approximation  to  the  test  probabilities. 
It  is  asymptotically  correct  and  performed  well  in  a  small  sample  simulation  study. 
The  test  is  based  on  the  maximum  norm  of  the  score  process  obtained  from  Co.x' 
partial  likelihood.  The  test  probabilities  are  obtained  by  a  ’parametric  bootstrap", 
i.e.  by  simulation  from  the  asymptotic  distribution,  with  an  unknown  variance  func¬ 
tion  replaced  by  an  estimate.  The  method  is  computationally  demanding,  but  still 
within  the  capabilities  of  a  standard  personal  computer.  An  important  advantage  is 
flexibility;  by  obvious  simple  changes  the  program  can  be  used  with  any  test  statistic 
based  on  the  score  process.  Some  problems  related  to  the  size  effect  in  the  strength 
of  materials  are  discussed,  and  the  method  is  applied  to  a  data  set  on  the  strengths 
of  carbon  fibers.  It  is  also  illustrated  on  two  cancer  studies  considered  by  Wei. 
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DONATAS  SURGAILIS 


Dr.  Donatas  Surgailis  of  the  Institute  of  Mathematics  and  Informatics  of  the 
Lithuanian  Academy  of  Sciences  in  Vilnius  visited  the  Center  for  three  months.  He 
introduced  a  new  rich  class  of  stationary  stable  processes  generalizing  moving  averages 
jointly  with  S.  Cambanis,  V.  Mandrekar  and  J.  Rosinski. 

1.  Generalized  stable  moving  averages  [1] 

No  explicit  representation  is  known  for  all  stationary  non-Gaussian  stable  pro¬ 
cesses.  The  main  two  subclasses  studied,  which  have  explicit  representations  moti¬ 
vated  by  the  Gaussian  case,  are  the  harmonizable  processes,  which  are  superpositions 
of  harmonics  with  stable  amplitudes,  and  the  moving  average  processes,  which  are 
filtered  white  stable  noise.  While  in  the  Gaussian  case,  the  latter  is  a  subclass  of 
the  former,  in  the  non-Gaussian  stable  case  the  two  classes  are  disjoint.  The  study 
of  stable  moving  average  processes  is  facilitated  by  the  fact  that  their  distribution 
is  essentially  (except  for  a  translation  and  sign)  determined  by  the  filter  function  of 
the  moving  average.  This  has  made  it  possible  to  study  distributional  properties  of 
the  process  (mixing,  ergodicity,  self-similarity,  Markov  property,  etc.)  through  the 
properties  of  the  filter  functions. 

In  this  work  the  class  of  non-Gaussian  stable  moving  average  processes  is  expanded 
substantially  by  the  introduction  of  an  appropriate  joint  randomization  of  the  filter 
function  and  of  the  stable  noise,  leading  to  stable  generalized  moving  averages  (GMA). 
The  characterization  of  their  distribution  through  their  filter  function  and  their  mixing 
measure  leads  to  a  far  reaching  generalization  of  a  theorem  of  Kanter  (1972). 

It  is  shown  that  stable  GMA’s  contain  sums  of  independent  stable  moving  averages 
and  that  they  are  still  disjoint  from  the  harmonizable  processes,  but  are  closed  under 
time  invariant  filters,  and  that  they  are  mixing,  so  they  have  strong  ergodic  properties. 
They  lead  to  a  wealth  of  new  examples  of  self-similar  processes,  beyond  the  linear 
fractional  stable  motions,  and  also  of  processes  which  are  reflection  positive,  which  is 
a  useful  weakening  of  the  Markov  property. 
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WEI  WU 


Professor  Wei  Wu  of  the  Statistics  Department  of  the  University  of  Illinois  visited 
the  Center  for  one  month  and  worked  jointly  with  S.  Cambanis  and  E.  Carlstein  on 
an  extensive  revision  and  generalization  of  a  part  of  her  Ph.D.  dissertation  research 
contained  in  the  following  paper. 


1.  Bootstrapping  the  sample  mean  for  data  with  infinite  variance  [1] 

When  data  comes  from  a  distribution  belonging  to  the  domain  of  attraction  of  a 
stable  law,  Athreya  (1987)  showed  that  the  bootstrapped  sample  mean  has  a  random 
limiting  distribution,  implying  that  the  naive  bootstrap  could  fail  in  the  heavy-tailed 
case.  The  goal  here  is  to  classify  all  possible  limiting  distributions  of  the  bootstrapped 
sample  mean  when  the  sample  comes  from  a  distribution  with  infinite  variance,  allow¬ 
ing  the  broadest  possible  setting  for  the  (nonrandom)  scaling,  the  resample  size,  and 
the  mode  of  convergence  (in  law).  The  limiting  distributions  turn  out  to  be  infinitely 
divisible  with  possibly  random  Levy  measure,  depending  on  the  resample  size.  An 
averaged-bootstrap  algorithm  is  then  introduced  which  eliminates  any  randomness 
in  the  limiting  distribution.  Finally,  it  is  shown  that  (on  the  average)  the  limiting 
distribution  of  the  bootstrapped  sample  mean  is  stable  if  and  only  if  the  sample  is 
taken  from  a  distribution  in  the  domain  of  (partial)  attraction  of  a  stable  law. 
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PH.D.  STUDENTS 


PH.D.  DEGREES  AWARDED 


DAVID  G.  BALDWIN 

Dr.  Baldwin  completed  his  Ph.D.  degree  working  under  the  direction  of  G. 
Kallianpur.  His  thesis  is  described  in  item  1  below. 

1.  Topics  in  the  theory  of  stochastic  processes  taking  values  in  the  dual 
of  a  countably  Hilbertian  nuclear  space  [1] 

•A.  theorem  is  given  on  the  weak  approximation  of  solutions  to  infinite  dimensional 
stochastic  differential  equations.  An  example  is  given  of  the  weak  approximation  of 
a  spatial  neuronal  model  with  reversal  potentials  by  a  continuous  diffusion  taking 
values  in  the  dual  of  a  countaly  Hilbertian  nuclear  space. 

Lastly  we  give  conditons  for  existence  and  uniqueness  of  global  McKean- Vdasov 
equations.  Results  are  extended  to  local  McKean-Vlasov  equations. 
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J.  XIONG 


Dr.  Xiong  completed  his  Ph.D.  degree  working  under  the  direction  of  G.  Kallian- 
pur.  His  thesis  is  described  in  item  1  below.  He  has  worked  jointly  with  G.  Kallianpur 
in  developing  various  aspects  of  the  theory  of  infinite  dimensional  stochastic  differen¬ 
tial  equations  listed  in  items  2  to  5  below. 

1.  Nuclear  space  valued  stochastic  differential  equations  driven  by  Poisson 
random  measures  [1] 

This  thesis  is  devoted  primarily  to  the  study  of  stochastic  differential  equations 
on  duals  of  nuclear  spaces  driven  by  Poisson  random  measures.  The  existence  of  a 
weak  solution  is  obtained  by  the  Galerkin  method  and  the  uniqueness  is  established 
by  implementing  the  Yamada-VVatanabe  argument  in  the  present  setup. 

When  the  magnitudes  of  the  driving  terms  are  small  enough  and  the  Poisson 
streams  occur  frequently  enough,  it  is  proved  that  the  stochastic  differential  equations 
mentioned  above  can  be  approximated  by  diffusion  equations. 

Finally,  we  consider  a  system  of  interacting  stochastic  differential  equations  driven 
by  Poisson  random  measures.  Let  (X”(f), . . . ,  X”(t))  be  the  solution  of  this  system 
and  consider  the  empirical  measures 

It  is  proved  that  converges  in  distribution  to  a  non-random  measure  which  is  the 
unique  solution  of  a  McKean-Vlasov  equation. 

The  above  problems  are  motivated  by  applications  in  neurophysiology,  in  partic¬ 
ular  to  the  fluctuation  of  voltage  potentials  of  spatially  distributed  neurons  and  to 
the  study  of  asymptotic  behavior  of  large  systems  of  interacting  neurons. 

2.  The  existence  and  uniqueness  of  the  solution  of  nuclear  space-valued 
stochastic  differential  equations  driven  by  Poisson  random  measures  (with 
G.  Hardy,  S.  Ramasubramanian  and  J.  Xiong)  [2] 

In  this  paper,  we  study  SDK’s  on  duals  of  nuclear  spaces  driven  by  Poisson  random 
measures.  The  existence  of  a  weak  solution  is  obtained  by  the  Galerkin  method. 
For  uniqueness,  a  class  of  £^-valued  processes  which  are  called  Good  processes  is 
introduced.  An  equivalence  relation  is  established  between  SDK’s  driven  by  Poisson 
random  measures  and  those  by  Good  processes.  The  uniqueness  is  established  by 
extending  the  Yamada-Watanabe  argument  to  the  SDK’s  driven  by  Good  processes. 
This  is  an  extension  to  discontinuous  infinite  dimensional  SDK’s  of  work  done  by  G. 
Kallianpur,  I.  Mitoma  and  P...  Wolpert  for  nuclear  space  valued  diffusions  [Stochastics. 
29,  1-45,  (1990)]. 

3.  Stochastic  differential  equations  in  infinite  dimensions:  A  brief  survey 
and  some  new  directions  of  research  [3] 

This  is  a  brief  survey  of  some  recent  work  on  nuclear  space  valued  stochastic 
differential  equations.  The  emphasis  is  on  stochastic  differential  equations  driven 
by  Poisson  random  measures.  An  application  of  the  evolution  equation  is  made  to 
stochastic  models  of  environmental  pollution.  The  asymptotic  behavior  of  interacting 
systems  of  nuclear  space  valued,  Poisson-driven  SDK’s  is  examined  and  a  propagation 
of  chaos  result  is  presented.  Some  new  directions  of  work  are  suggested. 
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4,  Asymptotic  behavior  of  a  system  of  interacting  stochastic  differential 
equations  driven  by  Poisson  random  measures  [4] 

In  this  paper,  we  study  a  system  of  interacting  stochastic  differential  equations 
taking  values  in  nuclear  spaces  and  driven  by  Poisson  random  measures.  VVe  also  con¬ 
sider  the  McKean- Vlasov  equation  associated  with  the  system.  VVe  show  that  under 
suitable  conditions  the  system  has  a  unique  solution  and  the  sequence  of  its  empirical 
distributions  converges  to  the  solution  of  the  McKean- Vlasov  equation  when  the  size 
of  the  system  tends  to  infinity.  The  results  are  applied  to  the  voltage  potentials  of 
a  large  system  of  neurons  and  a  law  of  large  numbers  for  the  empirical  measure  is 
obtained. 


5.  Stochastic  models  of  environmental  pollution  [5] 

In  this  paper,  we  consider  several  stochastic  models  arising  from  environmental 
problems.  First,  we  study  the  pollution  in  a  domain  where  undesired  chemicals  are 
deposited  at  random  times  and  locations  according  to  Poisson  streams.  Incorporated 
with  drift  and  dispersion,  the  chemical  concentration  can  be  modeled  by  a  linear 
stochastic  partial  differential  equation  (SPDE)  which  is  solved  by  applying  a  general 
result.  Various  properties,  especially  the  limit  behavior  of  the  pollution  process,  are 
discussed.  Secondly,  we  consider  the  pollution  problem  when  a  tolerance  level  is 
imposed.  The  chemical  concentration  can  still  be  modeled  by  a  SPDE  but  is  no 
longer  linear.  Its  properties  are  investigated  in  this  paper.  Finally,  the  linear  filtering 
is  considered  based  on  the  data  of  several  observation  stations. 
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DISSERTATION  IN  PROGRESS 


A.  BUDHIRAJA 

A  contribution  to  the  theory  of  McShane  stochastic  integrals 

Under  the  direction  of  G.  Kallianpur,  Amarjit  Budhiraja  is  developing  the  theory 
of  McShane  stochastic  integrals.  The  topic  contains  results  that  are  extensions  of 
some  of  the  results  obtained  by  G.W.  Johnson  and  G.  Kallianpur  and  reported  in  an 
earlier  annual  scientific  report. 
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[354]  S.  Cambanis,  Random  filters  which  preserve  the  normality  of  non-stationary 
random  inputs.  .Nov.  91.  Nonstationary  Stochastic  Processes  and  their  Applica¬ 
tions.  A.G.  .Miamee  ed..  World  Scientific.  1992,  219-237. 

[355]  J.  Olsson  and  H.  Rootzen,  .An  image  model  for  quantal  response  analysis  in 
perimetry,  .Nov.  91. 

[356]  Y.-C.  Su.  Sampling  designs  for  estimation  of  regression  coefficients  and  of  a 
random  process,  Dec.  91.  Dissertation. 

[357]  G.  Kallianpur,  Stochastic  differential  equation  models  for  spatially  distributed 
neurons  and  propagation  of  chaos  for  interacting  systems.  Dec.  91.  J.  .Math. 
Biol.,  to  appear. 

[358]  H.  Hurd  and  G.  Kallianpur,  Periodically  correlated  and  periodically  unitary 
processes  and  their  relationship  to  L^[0.  r]-valued  stationary  sequences.  Feb.  92. 
Nonstationary  Stochastic  Processes  and  their  .Applications.  .A.G.  .Miamee  ed.. 
World  Scientific.  1992,  256-287. 
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[359]  C.  Houdre,  Path  reconstruction  of  processes  from  missing  and  irregular  samples. 
Feb.  92.  Ann.  Probability,  to  appear. 

[360]  J.  Farshidi,  Autoregressive  expansion  of  the  linear  predictor  for  stationary 
stochastic  processes,  Mar.  92. 

[361]  D.  Monrad  and  H.  Rootzen,  Small  values  of  fractional  Brownian  motion  and 
locally  nondeterministic  Gaussian  processes.  Mar.  92. 

[362]  C.A.  Leon  and  J-C  Masse,  La  mediane  simpliciale  d’Oja:  existence,  unicite  et 
stabilite  ,  Mar.  92. 

[363]  1.  Fakhre-Zakeri  and  J.  Farshidi,  .A  central  limit  theorem  with  random  indices 
for  stationary  linear  processes,  Apr.  92. 

[364]  J.  Xiong,  Nuclear  space- valued  stochastic  differential  equations  driven  by  Pois¬ 
son  random  measures,  Apr.  92.  Dissertation 

[365]  D.  Surgailis,  J.  Rosinski,  V.  Mandrekar  and  S.  Cambanis,  Stable  generalized 
moving  averages,  Apr.  92. 

[366]  T.  Hsing,  Limit  theorems  for  stable  processes  with  application  to  spectral  den¬ 
sity  estimation,  June  92. 

[367]  G.  Kallianpur  and  V.G.  Papanicolaou,  Integration  over  Hilbert  spaces:  exam¬ 
ples  inspired  by  the  harmonic  oscillator,  July  92. 

[368]  H.  Hurd  and  A.  Russek,  Stepanov  almost  periodically  correlated  and  almost 
periodically  unitary  processes,  July  92. 

[369]  H.  Hurd  and  A.  Russek,  .Almost  periodically  correlated  processes  on  LC.A 
groups,  July  92. 

[370]  R.L.  Karandikar  and  V.G.  Kulkarni,  Second-order  fluid  flow  model  of  a  data- 
buffer  in  random  environment,  July  92. 

[371]  R.  Cheng,  Outer  factorization  of  operator  valued  weight  functions  on  the  torus, 
July  92. 

[372]  G.  Kallianpur  and  J.  Xiong,  Stochastic  differential  equations  in  infinite  dimen¬ 
sions:  A  brief  survey  and  some  new  directions.  Sept.  92. 

[373]  .A.G.  Bhatt,  G.  Kallianpur  and  R.L.  Karandikar.  On  interacting  systems  of 
Hilbert  space  valued  diffusions.  Sept.  92. 

[374]  C.  Houdre  and  A.  Kagan,  Variance  inequalities  for  functions  of  Gaussian  vari¬ 
ables,  Oct.  92. 

[375]  M.  Hernandez  and  C.  Houdre.  Disjointness  results  for  some  classes  of  stable 
processes,  Oct.  92. 

[376]  C.  Houdre,  Wavelets,  probability  and  statistics:  some  bridges.  Oct.  92. 

[377]  S.  Nandagopalan,  M.R.  Leadbetter  and  J.  Hiisler,  Limit  theorems  for  nonsta¬ 
tionary  strongly  mixing  vector  random  measures.  .Nov.  92. 
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[378]  R.  Lund,  A  dam  with  setisonal  input,  Nov.  92. 

[379]  R.  Cheng,  Operator  valued  functions  of  several  variables:  Factorization  and 
invariant  subspaces,  Nov.  92. 

[380]  I.  Fakhre-Zakeri  and  J.  Farshidi,  Limit  theorems  for  sample  covariances  of  sta¬ 
tionary  linear  processes  with  applications  to  sequential  estimation.  Sept.  92. 

[381]  J.  Farshidi,  G.  Kallianpur  and  V.  Mandrekar,  Spectral  characterization  and  au¬ 
toregressive  expansion  of  linear  predictors  for  second  ordered  stationary  random 
fields.  Part  1:  half  planes,  Dec.  92. 

[382]  J.  Farshidi,  Spectral  characterization  and  prediction  of  Z^’-representable 
stochastic  processes,  and  some  related  extremal  problems  in  ^'’-spaces,  (0  < 
p  <  oo),  Dec.  92. 

[383]  S.  Cambanis  and  I.  Fakhre-Zakeri,  On  prediction  of  heavy-tailed  autoregressive 
sequences:  Regression  versus  best  linear  prediction,  Dec.  92. 

[384]  A.  Deis  and  H.  Rootzen.  A  A:-sample  test  for  proportional  hazards  with  an 
application  to  the  strength  of  materials,  Dec.  92. 

[385]  H.  Rootzen,  Quantile  estimation  in  a  nonparametric  component  of  variance 
framework  with  applications  to  vision  problems,  Dec.  92. 

[386]  D.  Monrad,  Some  uniform  dimension  theorems  for  the  sample  functions  of  Levy 
processes  with  local  times,  Jan.  93. 


IN  PREPARATION 

K.  Benheni  and  S.  Cambanis,  The  effect  of  quantization  on  the  performance  of 
sampling  designs,  in  preparation. 

K.  Benhenni,  S.  Cambanis  and  Y.C.  Su,  Sampling  designs  for  regression  coeffi¬ 
cient  estimation  with  correlated  errors,  in  preparation. 

S.  Cambanis  and  C.  Houdre,  Wavelet  transforms  of  random  processes,  in  prepa¬ 
ration. 

G.  Kallianpur  and  R.L.  Karandikar,  Nonlinear  transformations  of  the  canonical 
Gauss  measure  on  Hilbert  space  and  absolute  continuity,  in  preparation. 

G.  Kallianpur  and  J.  Xiong,  .Asymptotic  behavior  of  a  system  of  interacting 
stochastic  differential  equations  driven  by  Poisson  random  measures,  in  prepa¬ 
ration. 

G.  Kallianpur  and  J.  Xiong,  Stochastic  models  of  environmental  pollution,  in 
preparation. 

M.R.  Leadbetter  and  H.  Rootzen,  Tail  estimation  for  stationary  sequences,  in 
preparation. 
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Seminars 


Sept.  11 
Sept.  19 
Sept.  25 

Oct.  2 

Oct.  9 

Oct.  11 

Oct.  16 

Oct.  18 

Oct.  23 

Oct.  30 

Xov.  6 
Nov.  13 

Nov.  20 

Nov.  21 

Nov.  26 


J.  StashefF,  UNC  -  Chapel  Hill:  Taking  the  measure  of  string  field 
theory 


Vasilis  G.  Papanicolaou,  Duke  University:  Multidimensional 
Schrodinger  operators  with  almost  periodic  potentials 


D.  Kolzow,  University  of  Erlangen- Niirnb erg:  Integral  transform  as¬ 
sociated  with  fractional  Brownian  motion 

C.  Houdre,  University  of  Maryland:  Path  reconstruction  of  processes 
from  irregular  samples 

K.  Petersen,  UNC  -  Chapel  Hill:  Random  sampling  of  stationary 
processes.  Including  a  disucssuion  of  the  divergent  views  of  Birkhoff, 
von  Neumann,  and  Wiener  on  convergence 


R.E.  Kalman,  ETH,  Zurich  and  University  of  Florida:  Identification 
of  linear  relations  in  noise 

V.  VVihstutz,  UNC  -  Charlotte:  Stochastic  averaging  and  stabilization 
by  random  vibration 

Y.  Kifer,  Hebrew  University,  Jerusalem:  Averaging  in  dynamical 
systems  and  large  deviations 

B.  Rozovskii,  University  of  Southern  California:  Stochastic  partial 
differential  equations  and  intermittency 

D.  Monrad.  University  of  Illinois  and  Center  for  Stochastic  Processes: 
Some  uniform  dimension  theorems  for  the  sample  functions  of  stable 
Levy  processes  with  local  times 

G.  Lawler,  Duke  University:  Nearest  neighbor  cluster  models 

R.  Andersen,  UNC  -  Charlotte:  Long  run  average  cost  for  an  optimal 
replacement  model  and  an  optimal  maintenance-replacement  model 
in  reliability 


J.-C.  Masse.  Laval  University  and  Center  for  Stochastic  Processes: 
On  Gaussian  reciprocal  processes 


T.  Norberg,  University  of  Gdteborg  and  Center  for  Stochastic 
Processes:  On  weak  lumpability  for  discrete  time  finite  state  space 
Markov  chains 

H.  Hurd:  Periodically  and  almost  periodically  correlated  random 
processes 
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Dec.  4 

Jan.  22 
Feb.  5 
Feb.  12 

Feb.  19 

Feb.  26 

Mar.  11 

Mar.  16 

Mar.  25 

.A.pr.  2 
Apr.  9 
-Apr.  15 
.Apr.  29 
Apr. -May 
May  1.3 
May  15 
May  20 


J.  Rissanen.  IBM  Research  Center,  San  Jose,  CA:  "Universal”  mod¬ 
eling  and  prediction  of  time  series 

J.A.  Cima,  UNC  -  Chapel  Hill;  Wavelets:  .4  theorem  of  Mallat 
J.A.  Cima,  UNC  -  Chapel  Hill:  Wavelets:  A  theorem  of  Daubechies 


S.I.  Resnick,  Cornell  University:  Convergence  of  scaled  random  sam¬ 
ples  in 

H.  Rootzen,  University  of  Lund  and  Center  for  Stochastic  Processes: 
What  is  the  probability  that  a  Gaussian  path  is  Bat? 

D.  Surgailis,  Lithuanian  Academy  of  Sciences,  Vilnius,  and  Center 
for  Stochastic  Processes;  Asymptotics  of  random  solutions  of  the 
Burgers  equation 

I.  Fakhre-Zakeri,  University  of  Maryland  and  Center  for  Stochastic 
Processes:  Models  of  Empirical- Bayes  type  for  software  reliability: 
IdentiBability  and  applications  to  optimal  stopping  of  software  testing 


B.V.  Rao.  Indian  Statistical  Institute  and  Indiana  University:  The 
dynamics  of  quadratic  maps  under  random  iteration 

A.  Russek,  Polish  Academy  of  Science,  Sopot:  Regularity  properties 
of  the  conditional  expectation  in  nonlinear  white  noise  Bitering 

I. M.  Sonin,  UNC-Charlotte:  The  asymptotic  behavior  of  Bnite  non- 
homogeneous  Markov  chains 

J. -.A.  Yan.  Institute  of  Applied  Mathematics.  Beijing:  White  noise 
calculus  in  terms  of  formal  series  expansions 

M.  Scarsini,  Universita  D’Annunzio.  Pescara.  Italy  and  Duke  Univer¬ 
sity;  Some  results  on  the  comparison  of  stochastic  processe 

P.  Hitczenko.  North  Carolina  State  University:  On  domination  in¬ 
equality  for  certain  martingale  transforms 

R.L.  Karandikar,  ISI  and  Center  for  Stochastic  Processes;  .4  short 
course  on  Markov  processes 

M.G.  Nadkarni.  University  of  Bombay  and  McGill  University:  Some 
spectral  questions  in  ergodic  theory 

R.M.  Gray.  Stanford  University:  Image  compression  and  vector  quan¬ 
tization:  Clustering  and  classiBcation  trees 


R.L.  Karandikar,  Indian  Statistical  Institute  and  Center  for  Stochas- 
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tic  Processes;  Invariant  measures  and  evolution  equations  for  Markov 
processes  characterized  via  martingale  problems 

May  27  T.  Zak,  Technical  University,  Wroclaw.  Poland:  Some  new  results  on 

Gaussian  measures;  Strictness  of  Anderson's  inequality  and  isoperi- 
metric  properties  of  strips 

June  3  M.  Burnashev,  Russian  Academy  of  Sciences,  Moscow:  On  a  search 

and  related  large  deviations  problems 

June  5  A.H.  Korezlioglu,  Institute  of  Telecommunications,  Paris:  Product 

form  approximations  of  finite  capacity  networks 

June  10  R.  Chaganty,  Old  Dominion  University,  Norfolk;  Large  deviations 

for  the  bootstrap  distrubutions 

June  17  I.  Fakhre-Zakeri,  University  of  Maryland  and  Center  for  Stochastic 

Processes:  .4  central  limit  theorem  with  random  indices  for  stationary 
linear  processes  with  applications 

June  24  C.  Houdre,  University  of  Maryland  and  Center  for  Stochastic  Pro¬ 

cesses:  Variance  inequalities  for  functions  of  Gaussian  variables 

July  1  R.L.  Karandikar,  Indian  Statistical  Institute,  Delhi,  and  Center  for 

Stochastic  Processes:  Stochastic  differential  equations  with  values  in 
a  Hilbert  space  and  propogation  of  chaos 

July  8  J.  Leskow,  University  of  California,  Santa  Barbara:  Inference  for 

repeatable  events 

July  15  J.  Farshidi,  Center  for  Stochastic  Processes:  Spectral  characteriza¬ 

tion  and  autoregressive  expansion  of  horizontal  and  vertical  linear 
predictions  for  stationary  second  order  random  fields 

July  22  R.  Cheng,  University  of  Louisville  and  Center  for  Stochastic  Pro¬ 

cesses:  A  Wold-type  decomposition  for  second-order  stationary  ran¬ 
dom  fields 

Aug.  6  I.  Weissman.  Technion,  Israel:  The  indices  of  the  largest  observations 

among  n  independent  ones 


PROFESSIONAL  PERSONNEL 

Faculty  Investigators 
S.  Cambanis 
G.  Kallianpur 
M.R.  Leadbetter 

Visitors 


R.  Cheng,  University  of  Louisville 

I.  Fakhre-Zakeri,  University  of  Maryland 

J.  Farshidi,  Michigan  State  University 
L.  Holst,  University  of  Stockholm 

C.  Houdre,  University  of  Maryland 
H.  Hurd 

R.L.  Karandikar,  Indian  Statistical  Institute, 
J.-C.  Masse,  University  of  Laval 

D.  Monrad,  University  of  Illinois 

T.  Norberg,  University  of  Goteborg,  Sweden 
J.  Olsson,  Lund  University,  Sweden 

V.  Papanicolaou,  Duke  University 
R.  Perfekt,  Lund  University,  Sweden 
H.  Rootzen.  Lund  University,  Sweden 

A.  Russek,  Polish  Academy  of  Science,  Sopot 
D.  Surgailis,  Lithuanian  Academy  of  Science, 

W.  Wu,  University  of  Illinois 


June  -  July  92 
January  -  December  92 
October  91  -  December  92 
January  92 
15  May  -  15  July  92 
entire  period 

New  Delhi  April  -  July  92 

mid  October  -  November  92 
September  -  December  91 
November  91 
November  91 
May  -  June  92 
September  -  October  91 
September  91  -  July  92 
January  -  December  92 
Vilnius  Sept.  -  Oct.  91,  February  92 

.August  92 


Graduate  Students 
D.  Baldwin 
A.  Budhiraja 
J.  Xiong 

Supported  by: 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 
NATIONAL  SCIENCE  FOUNDATION 

Faculty  Investigator:  S.  Cambanis 

Visitors:  R.  Cheng,  I.  Fakhre-Zakeri,  J.  Farshidi,  C.  Houdre,  D.  Monrad,  T.  Nor¬ 
berg,  H.  Rootzen,  D.  Surgailis,  VV.  Wu. 

Graduate  Student:  J.  Xiong 

ARMY  RESEARCH  OFFICE 

Faculty  Investigator:  G.  Kallianpur 


Visitors:  R.L.  Karandikar,  V.  Papanicolaou,  D.  Surgailis 


Graduate  Student:  J.  Xiong 
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